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Abstract 

We present a study of phi-four theory on noncommutative spaces using a combination 
of the Wilson renormalization group recursion formula and the solution to the zero di¬ 
mensional vector/matrix models at large N. Three fixed points are identified. The matrix 
model 9 = oo fixed point which describes the disordered-to-non-uniform-ordered transi¬ 
tion. The Wilson-Fisher fixed point at 9 = 0 which describes the disordered-to-uniform- 
ordered transition, and a noncommutative Wilson-Fisher fixed point at a maximum value 
of 9 which is associated with the transition between non-uniform-order and uniform-order 
phases. 

1 Introduction 

A noncommutative field theory is a non-local field theory in which we replace the or¬ 
dinary local point-wise multiplication of fields with the non-local Moyal-Weyl star prod¬ 
uct m- This product is intimately related to coherent states ms], Berezin quanti¬ 
zation [9] and deformation quantization m- It is also very well understood that the 
underlying operator/matrix structure of the theory, exhibited by the Weyl map [5], is 
the singular most important difference with commutative field theory since it is at the 
root cause of profound physical differences between the two theories. We suggest [3] and 
references therein for elementary and illuminating discussion of the Moyal-Weyl product 
and other star products and their relations to the Weyl map and coherent states. 

Noncommutative field theory is believed to be of importance to physics beyond the 
standard model and the Hall effect [Mj and also to quantum gravity and string theory 

[351136]. 
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Noncommutative scalar field theories are the most simple, at least conceptually, quan¬ 
tum field theories on noncommutative spaces. Some of the novel quantum properties of 
noncommutative scalar field theory and scalar phi-four theory are as follows: 

1. The planar diagrams in a noncommutative are essentially identical to the planar 
diagrams in the commutative theory as shown originally in [n]. 

2. As it turns out, even the free noncommutative scalar field is drastically different 
from its commutative counterpart contrary to widespread believe. For example, 
it was shown in [H] that the eigenvalues distribution of a free scalar field on a 
noncommutative space with an arbitrary kinetic term is given by a Wigner semicircle 
law. This is due to the dominance of planar diagrams which reduce the number of 
independent contractions contributing to the expectation value < 4>^'^ > from 2"'n! 
to the number iVpianar(2n) of planar contractions of a vertex with 2n legs. See 
also |45H48] for an alternative derivation. 

3. More interestingly, it was found in |12] that the renormalized one-loop action of a 
noncommutative suffers from an infrared divergence which is obtained when we 
send either the external momentum or the non-commutativity to zero. This non- 
analyticity at small momenta or small non-commutativity (IR) which is due to the 
high energy modes (UV) in virtual loops is termed the UV-IR mixing. 

4. We can control the UV-IR mixing found in noncommutative (j)'^ by modifying the 
large distance behavior of the free propagator through adding a harmonic oscillator 
potential to the kinetic term m- More precisely, the UV-IR mixing of the theory 
is implemented precisely in terms of a certain duality symmetry of the new action 
which connects momenta and positions [20]. The corresponding Wilson-Polchinski 
renormalization group equation |23p 24| of the theory can then be solved in terms of 
ribbon graphs drawn on Riemann surfaces. Renormalization of noncommutative 
along these lines was studied for example in HMS]. Other approaches to renormal¬ 
ization of quantum noncommutative (j)'^ can be found for example in |28H33j . 

5. In two-dimensions the existence of a regular solution of the Wilson-Polchinski equa¬ 
tion |23| together with the fact that we can scale to zero the coefficient of the har¬ 
monic oscillator potential in two dimensions leads to the conclusion that the standard 
non-commutative cj)'^ in two dimensions is renormalizable m- In four dimensions, 
the harmonic oscillator term seems to be essential for the renormalizability of the 
theory m- 

6. The beta function of noncommutative (j)'^ theory at the self-dual point is zero to all 
orders |251127] . This means in particular that the theory is not asymptotically free in 
the UV since the RG flow of the coupling constant is bounded and thus the theory 
does not exhibit a Landau ghost, i.e. not trivial. In contrast the commutative (j)'^ 
theory although also asymptotically free exhibits a Landau ghost. 

7. Noncommutative scalar field theory can be non-perturbatively regularized using ei¬ 
ther fuzzy projective spaces CP"" |49j or fuzzy tori T” [43]. The fuzzy tori are inti¬ 
mately related to a lattice regularization whereas fuzzy projective spaces, and fuzzy 
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spaces [661167] in general, provide a symmetry-preserving sharp cutoff regularization. 
By using these regulators noncommutative scalar field theory on a maximally non¬ 
commuting space can be rewritten as a matrix model given by the sum of kinetic 
(Laplacian) and potential terms. The geometry in encoded in the Laplacian in the 
sense of [521E3]- 

The case of degenerate noncommutativity is special and leads to a matrix model 
only in the noncommuting directions. See for example [68] where it was also shown 
that renormalizability in this case is reached only by the addition of the doubletrace 
term f d^x{Tr4>)‘^ to the action. 

8. Another matrix regularization of non-commutative can be found in liEIlES] 
where some exact solutions of noncommutative scalar field theory in background 
magnetic fields are constructed explicitly. Furthermore, in order to obtain these 
exact solutions matrix model techniques were used extensively and to great efficiency. 
For a pedagogical introduction to matrix model theory see |69H74j . Exact solvability 
and non-triviality is discussed at great length in m- 

9. A more remarkable property of quantum noncommutative is the appearance of 
a new order in the theory termed the striped phase which was first computed in a 
one-loop self-consistent Hartree-Fock approximation in the seminal paper I3Z|. For 
alternative derivations of this order see for example [381139] . It is believed that the 
perturbative UV-IR mixing is only a manifestation of this more profound property. 
As it turns out, this order should be called more appropriately a non-uniform ordered 
phase in contrast with the usual uniform ordered phase of the Ising universality 
class and it is related to spontaneous breaking of translational invariance. It was 
numerically observed in d = 4 in [40| and in d = 3 in [mill] where the Moyal-Weyl 
space was non-perturbatively regularized by a noncommutative fuzzy torus [43] . The 
beautiful result of [41] shows explicitly that the minimum of the model shifts to a 
non-zero value of the momentum indicating a non-trivial condensation and hence 
spontaneous breaking of translational invariance. 

10. Therefore, noncommutative scalar 0“^ enjoys three stable phases: i) disordered (sym¬ 
metric, one-cut, disk) phase, ii) uniform ordered (Ising, broken, asymmetric one-cut) 
phase and hi) non-uniform ordered (matrix, stripe, two-cut, annulus) phase. This 
picture is expected to hold for noncommutative/fuzzy phi-four theory in any dimen¬ 
sion, and the three phases are all stable and are expected to meet at a triple point. 
The non-uniform ordered phase [54] is a full blown nonperturbative manifestation 
of the perturbative UV-IR mixing effect m which is due to the underlying highly 
non-local matrix degrees of freedom of the noncommutative scalar field. In [371138] . 
it is conjectured that the triple point is a Lifshitz point which is a multi-critical point 
at which a disordered, a homogeneous (uniform) ordered and a spatially modulated 
(non-uniform) ordered phases meet [55] . 

11. In [38| the triple (Lifshitz) point was derived using the Wilson renormalization group 
approach [56] . where it was also shown that the Wilson-Fisher fixed point of the 
theory at one-loop suffers from an instability at large non-commutativity. See [631164] 
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for a pedagogical introduction to the subject of the functional renormalization group. 
The Wilson renormalization group recursion formula was also used in |581162j to 
study matrix scalar models which, as it turns out, are of great relevance to the limit 
9 —> oo of noncommutative scalar field theory |65| . 

12. The phase structure of non-commutative in d = 2 and d = 3 using as a regulator 
the fuzzy sphere was studied extensively in |76H82j . It was confirmed that the phase 
diagram consists of three phases: a disordered phase, a uniform ordered phases 
and a non-uniform ordered phase which meet at a triple point. In this case it is 
well established that the transitions from the disordered phase to the non-uniform 
ordered phase and from the non-uniform ordered phase to the uniform ordered phase 
originate from the one-cut/two-cut transition in the quartic hermitian matrix model 
[691I7D]. The related problem of Monte Carlo simulation of noncommutative (j)^ on 
the fuzzy disc was considered in [83] . 

13. The above phase structure was also confirmed analytically by the multitrace ap¬ 
proach of [841185] which relies on a small kinetic term expansion instead of the usual 
perturbation theory in which a small interaction potential expansion is performed. 
This is very reminiscent of the Hopping parameter expansion on the lattice [87]I88]. 
See also [86] for a review and an extension of this method to the noncommutative 
Moyal-Weyl plane. For an earlier approach see [Sj and for a similar more non- 
perturbative approach see [IHtilH] . This technique is expected to capture the matrix 
transition between disordered and non-uniform ordered phases with arbitrarily in¬ 
creasing accuracy by including more and more terms in the expansion. Capturing 
the Ising transition, and as a consequence the stripe transition, is more subtle and 
is only possible if we include odd moments in the effective action and do not impose 
the symmetry (p —> —p. 

14. The multitrace approach in conjunction with the renormalization group approach 
and/or the Monte Carlo approach could be a very powerful tool in noncommutative 
scalar field theory. For example, multitrace matrix models are fully diagonalizable, 
i.e. they depend on N real eigenvalues only, and thus ergodic problems are absent 
and the phase structure can be probed quite directly. The phase boundaries, the 
triple point and the critical exponents can then be computed more easily and more 
efficiently. Furthermore, multitrace matrix models do not come with a Laplacian, 
yet one can attach to them an emergent geometry if the uniform ordered phase 
is sustained. See for example [96]. Also, it is quite obvious that these multitrace 
matrix models lend themselves quite naturally to the matrix renormalization group 
approach of [901193] . 

Among all the approaches discussed above, it is strongly believed that the renormal¬ 
ization group method is the only non-perturbative coherent framework in which we can 
fully understand renormalizability and critical behavior of noncommutative scalar field 
theory in complete analogy with the example of commutative quantum scalar held the¬ 
ory outlined in [89] . The Wilson recursion formula, in particular, is the oldest and most 
simple and intuitive renormalization group approach which although approximate agrees 
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very well with high temperature expansions [56]. In this approximation we perform the 
usnal truncation but also we perform a reduction to zero dimension which allows explicit 
calculation, or more precisely estimation, of Feynman diagrams. 

The goal in the first part of this article is to apply this method to scalar field 
theory at the self-dual point on a degenerate noncommutative spacetime with two strongly 
noncommnting directions. See also m- In the matrix basis this theory becomes, after 
appropriate non-perturbative definition, an N xN matrix model where is a regulator in 
the noncommutative directions, i.e. N here has direct connection with noncommutativity 
itself. More precisely, in order to solve the theory we propose to employ, following [58II60] , 
a combination of 

• i) the Wilson approximate renormalization group recursion formula 
and 

• ii) the solution to the zero dimensional large N counting problem given in this case 
by the Penner matrix model which can be tnrned into a mnltitrace matrix model for 
large values of 6. 

As discussed neatly in [58| the virtue and power of combining these two methods lies in 
the crucial fact that all leading Feynman diagrams in 1/N will be counted correctly in 
this scheme including the so-called ’’setting sun” diagrams. As it turns out the recursion 
formula can also be integrated explicitly in the large N limit which in itself is a very 
desirable property. 

In the second part of this article a non perturbative study of the Ising universality 
class fixed point in noncommutative 0{N) model is carried out using precisely a combi¬ 
nation of the above two methods. See also |94] . It is found that the Wilson-Fisher fixed 
point makes good sense only for sufficiently small values of 0 up to a certain maximal 
noncommutativity. This fixed point describes the transition from the disordered phase to 
the uniform ordered phase in the same way that the matrix model fixed point, obtained 
in the first model, describes the transition from the one-cut (disordered) phase to the 
two-cut (non-uniform ordered, stripe) phase. 

Another fixed point termed the noncommutative Wilson-Fisher fixed point is identified 
in this case. It interpolates between the commutative Wilson-Fisher fixed point of the 
Ising universality class which is found to lie at zero value of the critical coupling constant 
a* of the zero dimensional reduction of the theory and a novel strongly interacting fixed 
point which lies at infinite valne of a* corresponding to maximal noncommutativity. This 
is identified with the transition between non-uniform and uniform orders. 

This article is organized as follows: 

1. The 9 = oo Fixed Point in Self-Dual Degenerate Noncommutative 

• Degenerate Noncommutativity. 

• Wilson RG Recursion Formula. 

• The Zero-Dimensional Matrix Model. 

• IjO Correction. 

• 9 = oo Fixed Point and Critical Exponents. 
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• On the Wave Function Renormalization. 

2. The 0 = 0 Fixed Point in Noncommutative 0{N) Sigma Model. 

• Maximally Noncommuting 0{N) Sigma Model. 

• The Noncommutative Wilson-Fisher Fixed Point. 


2 The 9 = oo Fixed Point in Self-Dual Degenerate 
Noncommutative 

2.1 Degenerate Noncommutativity 

We are interested in phi-four theory on a degenerate noncommutative Moyal-Weyl 
space Rg = x R^ with a harmonic osicllator term is give by the action 


S 



-df + - df. H- 

2 2 2 ^2 




( 2 . 1 ) 


In terms of operators this reads 




S = 1^2 I dPx Tv'n 


-^ 9 ! + - ^52 + ^ ) $ + 5$^ 


( 2 . 2 ) 


The index i runs over the noncommuting directions while the index fi runs over the 
commuting directions. The Planck volume 1^2 is defined by 1^2 = 27r0 where 9 is the 
noncommutativity parameter and Xi = 2{6~^)ijXj. The parameters of the model are the 
mass rn^, the quartic coupling g, the harmonic oscillator parameter 9^. We can expand the 
scalar fields in the Landau basis {(pm,n} as (with x standing for commuting coordinates) 


6 = 


1 


OO 

m,n=l 


(2.3) 


Furthermore, by introducing a matrix regularization we obtain the action (with g = 
V2u/N) 


s = 


= J d^xTxM + a[EM‘^ + V^F+MFMJ 


(2.4) 


The coupling constants of the theory are the mass m?, the quartic coupling constant 
u/N, the noncommutativity parameter 9 and the harmonic oscillator parameter 9^. The 
parameters a and \/tJ are dehned by 


92 + 1 _ 92-1 


(2.5) 
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The external sources E and T are the matrices given by 


(r)zm — 15lm—l , {E)lm — {I 

At the self-dual point we have 17^ = 1 and thus the theory becomes 
5 = 


xTtn 


+ aEM'^ 


( 2 . 6 ) 


(2.7) 


2.2 Wilson RG Recursion Formula 

The Wilson renormalization group approach consists in general in the three main steps; 
1) Integration, 2) Rescaling and 3) Normalization. In our case here we will supplement the 
first step of integration with two approximations a) Truncation and b) Wilson Recursion 
formula. 


Integration: We start by decomposing the N xN matrix M into an N x background 
matrix B and an x fluctuation matrix E, viz M = B + E. The background B 
contains slow modes, i.e. modes with momenta less or equal than pA while the fluctuation 
E contains fast modes, i.e. modes with momenta larger than pA where 0 < p < 1. The 
integration step involves performing the path integral over the fluctuation E to obtain an 
effective path integral over the background B alone. We find 




j dB exp(-5[R] - A5(R)). 


( 2 . 8 ) 


An exact formula for AS{B) up to the fourth power in the field B is given by the cumulant 
expansion 


AS{B) 


+ 


4^ / < TtnB^F^x) >co 

8^ J d^x J d^y < TrNBF\x).TrNBF\y) >eo 

8 ^ / J d'^y < TrNB^F\x).TrNB^F\y) >,o 



d^z < TrNBF\x).TrNBF\y).TrNB^F\z) >co 

d^z [ d^w < TrNBF^{x).TrNBF^{y).TrNBF^{z).TrNBF^{w) >co ■ 


(2.9) 


The contribution in the ith line will be denoted ASi{B) in the following. The notation 
”co” stands for the connected component. The first and second terms yield correction to 
the mass parameter whereas the last three terms yield correction to the quartic coupling 
constant. The wave function renormalization is obtained from the expansion around 
= 0 of the second term which is the most difficult contribution to calculate. 
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The formula (j2.9ll is still very complicated. To simplify it and to get explicit equa¬ 
tions we employ the so-called Wilson truncation and Wilson recursion formula. This is 
usually thought of as part of the integration step. Wilson truncation means that we 
calculate quantum corrections to only those terms which appear in the original action. 
Wilson recursion formula is completely equivalent to the use in perturbation theory of the 
Polyakov-Wilson rules given by the following two rules: 


1) We replace every internal propagator 1/(A;^ -|- /U^) by 1/(A^ -|- m?) 
and 

2) We replace every momentum loop integral by another constant 

vd = where ud is given by 

2(11 1 


VD = 


D (47r)^/2r(Zl/2)' 


( 2 . 10 ) 


This is a very long and tedious calculation. The results, in the limit 9 —> oo, are as 
follows. 

Quantum corrections to the mass parameter and to the harmonic oscillator coupling 
constant a are obtained from the first term of (12.91) and also from the second term of (12.9p 
evaluated at kP' = 0. The corresponding Feynman diagrams are shown on figures ([T]) and 
([2]). We get 


'S'mass-l-harmonic oscillator — j d xTTj\[ 


+ aEB‘^ + ASi(5) + A52(5)|fc2=o 


d^xTrjsf -|- Attiq -|- oAm^ -|- 0{oP)')B‘^ -|- a(Aao -|- aAui -|- 0{aP)^EB‘^ 

( 2 . 11 ) 


The wave function renormalization is also obtained from the 2nd term of (12.9p and as a 
consequence the relevant Feynman diagrams are still given by those shown on figure ([2]). 

More precisely we need, as before, to expand these diagrams around k^ = 0 but retain 
now the linear term in k"^ which is very difficult to do explicitly. Our estimation of the 
coefficient of , motivated by dimensional consideration, is obtained by the approximation 
of [SHIEQ]. Explicitly we have the (3rd) rule: 

• 3) We approximate the first derivative of the propagator with respect to the external 
momentum k'^ by the multiplication with the given propagator as follows 

+ + ( 2 . 12 ) 

We get then 

Skinetic + dk2AS2y=o = [^ + Z + aAZ + 0{a^)] J d^xTr^id^Bf. (2.13) 

The renormalization of the quartic coupling is obtained from the last three terms of 
equation (|2.9p . i.e. from A^s, A54 and A^s, and is given explicitly 

‘S'interaction + A5'3|^2 _q -|- AS' 4 |;;. 2 _q -|- AS' 5 |^ 2 _g — jv Auq -I- uAui -|- 0(o^)) / d^xTr]\fB^[x). 








The quantum corrections AtUq, Amf, Auq, Aai, Z, AZ, Auq and Aui will be given 
explicitly in the next section. 

Scaling and Normalization: By performing the second step of the Wilson renor¬ 
malization group approach, i.e. by scaling momenta as p —)• p/p so that the cutoff returns 
to its original value A and the third and final step of the Wilson renormalization group 
approach consisting in rescaling the field in such a way that the kinetic term is brought 
to its canonical form we obtain the effective action 


m 


S + AS = - xTr]\f{d^B Y -\ -/ d^xTr^B^ + a / d^xTr^EB 


+ 


2 „ 
/ 
u 

N 


d^xTvNB'^. 


(2.15) 


The renormalized field B is related to the bare field B as follows. If B and B are the 
Fourier transforms of B and B respectively then 


B{p) = + 2(Z + aAZ) + 0{a?)13{pp). 


(2.16) 


The renormalized mass m ^, the renormalized quartic coupling constant u and the renor¬ 
malized inverse noncommutativity a are given by (with e = A — D) 


_ 2 m^ -|- Amg -|- aAm\ + 0{a‘^) 

(2.17) 

1+ 2(Z +aAZ)+ 0(a2) ' 

_2 a(Aao + aAai -|- 0(a^)) 

(2.18) 

l + 2(Z + aAZ) + 0(a2)' 

gtt(l -|- Attg -|- aAui + 0(a^)) 

(2.19) 

(1 + 2(Z + aAZ) + 0(a2))^’ 


The process which led from the bare coupling constants a and u to the renormalized 
coupling constants m'^, a"^ and u can be repeated an arbitrary number of times. The bare 
coupling constants will be denoted by mg, uq and uq whereas the the renormalized coupling 
constants at the first step of the renormalization group procedure will be denoted by mf, 
ai and ui. At a generic step Z -|- 1 of the renormalization group process the renormalized 
coupling constants m^^^, a;+i and are related to their previous values mf, ai and ui 
by precisely the above renormalization group equations. We are therefore interested in 
renormalization group flow in a 3—dimensional parameter space generated by the mass 
m?, the quartic coupling constant u and the harmonic oscillator coupling constant (inverse 
noncommutativity) a. 
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2.3 The Zero-Dimensional Matrix Model 

The explicit calculation of the corrections Atuq, Z and Auq, at 9 = oo using the above 
rules, reduces to the properties of the zero-dimensional matrix model 

r = + As*] , 9 = (2.20) 

The Schwinger-Dyson identity of this model can be rewritten in terms of the Green’s 
functions G 2 =< Tr^B'^ > /iV^ and G 4 =< Tr^B^ > /N^ as 

1 = G2 + 45G4. (2.21) 

The model ()2.20p is exactly solvable. The connected 2—point and 4—point functions C 2 
and C 4 and the 2—point and 4—point proper vertices r 2 and r 4 of this model are given 


by (with = (yTTdS^ — l)/2‘ig) 

C2 = G2,T2 = (C2)-^ , G 2 = - r2). (2.22) 

C4 = G4 - 2(G2)2 , r4 = -C4(C'2)-^ , Gi = r\3 - r^). (2.23) 

Thus, the functions r2(5') and r4(5() are known non-perturbatively given by 

^2(5) = ^2(41^2) = 1 + 85 - 8O5' + 16645^ - .... (2.24) 

r4(g) = =4g- 32ff^ + 896ff3 + ... (2.25) 

r^(4 — 


The Schwinger-Dyson identity of this model can also be rewritten in terms of r 2 and r 4 
as 


r2 — 1 -|- 8 gG2 — 4g'r4(G2)^. (2.26) 

The corrections Atuq and Auq are found to be given in terms of the 2—point proper vertex 
r 2 ( 5 ') and the 4—point vertex r 4 (gf) of the above zero-dimensional matrix model by 

Amg = (A^-h m^)(r2(5f) - 1). (2.27) 

Auo = - 1. (2.28) 

45 

Similarly, the wave function renormalization Z is found perturbatively to be given by the 
expansion 

z = 8g‘^ - 256^ -F ... (2.29) 

We need now to find a combination of Green’s functions and proper vertices of the above 
zero-dimensional matrix model with an expansion given exactly by 8 g'^ — 256g^ + .... 
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From the Schwinger-Dyson identity of the model (12.201) we propose that the function 
2gV/^{g)G 2 {g) is the correct guess. Notice the resemblance of the graphs corresponding to 
SgG 2 and —AgT^G^ to the graphs associated with the terms AS"! and AS '2 respectively. 
Indeed, we compute 

Z = 2gTMGl{g) 

2g r^(l — r^)(5 — 2r^) 

Y 4-r2 

= 85 ^- 2565 ^ + ... (2.30) 


The renormalization group equations are therefore given by 

'2 _2 m^ + 2 (m^ + A^)(r2(g) - 1 ) / _ ^A{g)u _ 

jjl — n _iii_— 'll = n _-_-_. 

1 + 4gT4{g)G2{g) 4^(1 + 4gT4{g)G2{g)Y 


(2.31) 


2.4 1/0 Correction 

We start with two remarks: 


1. The free propagator of this theory is simple given by 

/j .2 _j_ j • 


(2.32) 


In the limit 9 —> 00 this propagator behaves as 1/(A:^ + m^). In the computation 
of the effective action we need extensively the sum Yli For 9 = 00 this sum 

is obviously of order N. Including also the sub leading 1/9 correction this sum takes 
then the following form 


Y/ ^ij{k) —^ NAnojik) , no = ^ ^ ■ (2.33) 


A straightforward generalization of this result is 


i 

Again in the spirit of the Wilson contraction we will need to treat the index j in 
the propagator An^j (A) as a continuous variable and expand the propagator around 
j = n where n is some index. This actually makes sense since we are assuming that 
9 is sufficiently large and thus a is sufficiently small. Similarly to the expansion 
around = 0 , only the first two terms in the expansion around j = n are relevant 
to renormalization here. We choose n = 1/2 because the harmonic oscillator term 
is of the form f d^xB^(x)u(i — 1/2). From these considerations We have then the 
extra (4th) rule 

• 4) We expand An^j (A) around j = n as 

Anoj(c) = A„p„(A) — (j — n)aA^^„(A) + (j — n) a + ...(2.35) 
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2. The explicit calculation of the various quantum corrections using the above rules, 
for 9 7 ^ oo, reduces to the properties of the zero-dimensional matrix model 

V = TrN[\B^ + %B^ + A^EM^]. (2.36) 

2 JS 

This model we do not know how to solve exactly, similarly to the a = 0 model, and 
thus our results below will be given as perturbative expansions in g. 

The corrections Aiuq, Z and Auq are still given by the results of the previous section with 
the redefinition of g as 


UVD 

^ {A? + 171^ + aNy ' 


(2.37) 


On the other hand, the corrections Ami, Aao, Aoi, AZ and Am are given by the 
perturbative expansions 


Ami = N{r2{g) - Ar2{g)) , Ar2(<7) = 1 - 4^ + - 22405^ + .... (2.38) 

Aao = Ar2(5) , Aai = ^(^2 + ^ 2 ) - 80g^ + 512 / + ...). (2.39) 


AZ 


1024/+ 


(2.40) 


Aui 


N 

2 {K^ + m?) 


(-16(7 + 832/ + ...).. 


(2.41) 


2.5 9 = 00 Fixed Point and Critical Exponents 

By definition a renormalization group fixed point is a point in the space parameter 
which is invariant under the renormalization group flow. If we denote the fixed point by 
ml, a* and u* then we must have 


(m^ + Arug + aAmI + 0(a^))* 

(2.42) 

(l + 2(Z + aAZ) + 0(a2))* ’ 

2 a*(Aao + aAai + 0(a^))* 

(2.43) 

(l + 2(Z + aAZ) + 0(a2))/ 

m(l + Auo + aAui + 0(a2))* 

(2.44) 

(l + 2(Z + aAZ) + 0(a2))2 ' 


The second equation is new by comparison with the commutative theory. The definition 
of g^ in terms of ml, a* and u* is obvious. There are possibly several soultions (fixed 
points) of interest to these renormalization group equations. We will mainly concentrate 
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on the matrix model fixed point corresponding to infinite noncommutativity which is the 
most obvious solution to equation (j2.43j) given by 


a* = 0. 


The remaining two equations reduce then to 


2 -2 
m^= p 


{m? + Amp)* 

(1 + 2 ^)* 


(2.45) 


(2.46) 


U^: = P 


n*(l + Auo)* 

(1 + 2 Z )2 


(2.47) 


Thus this fixed point is fully determined by functions which are known non-perturbatively. 
An obvious solution to (I2.47P is u* = 0 which corresponds to the usual Gaussian fixed 
point. By discarding this solution equation (|2.47l) becomes 


^ ^ 4g,{l + 2Z{g,)y 

The critical value of the mass parameter is obtained from equation (I2.46p 

ml ^ p-‘^{T2{g*) - 1 ) 

A2 l + 2Z{g,)-p-^r2ig,)' 

The physical region of g is [0, oo[ while the full domain of definition is [- 


(2.48) 


(2.49) 
1/48, oo[. Fur- 


Did) 

g* 

ml/ 


1(3) 

7.603 

-0.935 

0.204 

2(4) 

2.282 

-0.851 

0.854 

3(5) 

0.409 

-0.643 

3.527 

4(6) 

0 

0 

0 


Table 1: The critical values for p = 1/2. 


thermore, the functions r 2 , r 4 and Z 2 = 3Z/2g depend on g only through r = r{g) defined 
by = 2 /(vT^F485 -|- 1). Graphically we observe that the above equation (I2.47P admits 
a non-trivial solution for all dimensions Zl = 1, 2,3,4 corresponding to d = 3,4,5, 6. The 
numerical solution for g^,, n* and ml are shown on table ©• There is of course in each 
dimension the extra Gaussian fixed point as we have discussed. There is only the Gaussian 
fixed point in D = 4(d = 6) in this approximation. Also, in our approximation we have 
checked that there is always a non-trivial fixed point for any value of p in the interval 

0 < p < 1. 

In the remainder we will compute the mass critical exponent and the anomalous 
dimension g within this scheme. 
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The computation of the mass critical exponent v, and other critical exponents, requires 
linearization of the renormalization above group equations. These renormalization group 
equations are of the form 

= M{G^^\p). (2.50) 


The vector of coupling constants G is defined by G = (Gi,G 2 ,G 3 ) where Gi = m?, 
G 2 = u and G 3 = a. The linearized renormalization group equations are of the form (with 
5G = G-G,) 

6 G^^+^^ = M{G^,p)6G^^\ (2.51) 


In our problem the matrix M is of the form 

^ Mil Mi2 Mi3 \ 

M21 M22 M23 1 

\ 0 0 M 33 j 


The eigenvalue in the direction G 3 = a is therefore given by 


A 3 — M 33 


p ^(Aap)^ 

1 + 2Z* 


(2.52) 


(2.53) 


This eigenvalue is plotted on figure ([3]) as a function of In p. It looks like that a is an 
irrelevant coupling constant. However, the function Aao used in the above formula is only 
known perturbatively and hence this conclusion should be taken with care. 

The two remaining eigenvalues are determined from the linearized renormalization 
group equations in the 2—dimensional space generated by Gi = and G 2 = u. These 
are given by 


r.].+% 


m 


dr 2 


2m^ dZ 
1 + 2Z dg 


(2.54) 


5u = 5u + 5gutf 


1 


d M 4 


(-)- 
^ A n ' 


dZ 


1 + Arto dg^Ag^ 1 + 2Z dg 


(2.55) 


6 g = —6u — „ (2.56) 

tt* A2 + m2 

As it turns out this problem depends only on functions which are fully known non- 
perturbatively. The eigenvalues Ai and A 2 can be determined from the trace and de¬ 
terminant which are given by 


Ai -|- A 2 — Mil + M 22 = Tr 2 M , A 1 A 2 — M 11 M 22 — M 12 M 21 = det 2 M. (2-57) 


In other words 


Ai — 


Tr2M ± v^(Tr2M)2 - 4det2M 


A 2 — Tr2Af — Ai 


(2.58) 


2 
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The eigenvalues Aj must scale with the dilatation parameter p as 


\i{p) = \i{\)p-y^. (2.59) 

The exponents pi are called critical exponents or scaling indices. The mass critical expo¬ 
nent V is given by the inverse of the critical exponent of the largest eigenvalue. If Ai > A 2 
then 


1 / = 1 / 2 / 1 . (2.60) 

We found that the solutions Ai and A 2 exist for 0.1 < /? < 1 for D = 3 and for 0.35 < p < 1 

for L) = 2 with the property Ai > A 2 > 0. The formula ()2.59l) was used then as a crucial 

test for our numerical calculations. In particular, we have determined by means of this 
formula the range of the dilatation parameter p over which the logarithm of the eigenvalues 
scale linearly with Inp. The eigenvalue InAi was found to be linear over the full range 
whereas the eigenvalue lnA 2 was linear only for Inp << — 1. In any case, we expect the 
behavior (I2.59P to hold only if the renormalization group steps are sufficiently small so 
not to alter drastically the infrared physics of the problem. Some results are shown on 
figure dH). We find explicitly for InAi the following fits: 

D = 2 , In Ai = -1.150(10) In p + 0.077(5) , — = 1.150(10). (2.61) 

D = 3, In Ai =-1.465(3) In p-F 0.027(2) , — = 1.465(3). (2.62) 

We conclude immediately that the scaling field corresponding to the mass is relevant since 
yi = l/vi > 0 while the scaling field corresponding to the quartic coupling constant is 
irrelevant, i.e. yi = 2 / 1^2 < 0, as seen immediately from the behavior of the eigenvalues 
for In p < < 1 on figure (jl]) . We skip writing explicitly the corresponding estimate of the 
coupling constant critical exponent 2 / 2 - 

In order to compute the anomalous dimension tj we need to go back to the wave 
function renormalization contained in equation ()2.16l) which can be put in the form 


-/ 2+D-r? ~ 

B{p) = p 2 B{pp). 

The coefficient p is called the anomalous dimension. It is given explicitly by 


(2.63) 


hl(l -|- 2 Z^:) 

V = -]- 

In p 

^ e _ ln(r4(ff^)/4ff^) 
2 2 Inp 


(2.64) 


The results in 74 = 3 and D = 2 are shown on figure ([5]). We observe that p approaches 
a constant value as Inp —> —1 for D = 2. 


15 











2.6 On the Wave Function Renormalization 


The wave function renormalization (|2.12l) can be improved by replacing the overall 
minus sign multiplying this equation by the correct coefficient coming from the leading 
Feynman diagrams. This coefficient is conjectured in [95] to be the same for all other 
subleading Feynman diagrams. The consequences of this change on the fixed point and 
the critical exponents can be found in |95] . 


3 The ^ = 0 Fixed Point in Noncommutative 0{N) 
Sigma Model 

3.1 Maximally Noncommuting 0{N) Sigma Model 

The action of interest in this section is of the form 

S = J (fx^ai-df + + Sint + U j (fx (^>0 * ^>a)^ + V j d'^X ($a * ^>6)^. (3.1) 

The vertex is given explicitly in momentum space by 

F(A:i,A: 2,A:3,A:4) = wcos-^l^cos^^l^ + |(e5(^i^^3+fc2Afc4) +g|(fciAfc4+fc2Afc3))^ (3 2) 

We decompose the fields into backgrounds (j)a{x) which contain slow modes, i.e. 

modes with momenta less or equal than pA and fluctuations faix) which contain fast 
modes, i.e. modes with momenta larger than pA where 0 < p < 1. The partition function 
is then given by 


Z = 






(3.3) 


By using the symmetry under (j)a —> —4>a and momentum conservation, we compute up 
to the 4th order in the slow fields cjta the non perturbative expansion 


^ AA 2 >co 2 ^ AAi >co 2 ^ -^^2 ^co T2 ^ AAiAA.2 >co ^ AA^ ^co • 

(3.4) 

At this stage we employ the large N limit. After appropriate rescaling, the propagator 
comes with 1/N factor, the vertex comes with a factor of N and the contraction of a vector 
index yields a factor of N. There exists a non-trivial 1/N expansion only if u,v —)• 0 
when N —)• 00 such that uq = uN and vo = vN is kept fixed. By inspection it is found 
that all terms of the form * fa* ■■ are subleading in the large N limit. In other 

words, we can set in the large N limit Aii,Ai 3 —)• 0 and 


M 2 = 2u dl’-x (t>a* 4>a* h* fb + 2v / dl’-x 4>a* fb*4>a* fb 
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(3.5) 






As a result the final form of the effective action is given explicitly by the very simple 
cumulant expansion 

A5eff[</>] = < M 2 >co -^ < Ml >co ■ (3.6) 

Next we will give the exact solution of the model in the large N limit by computing 
formally all Feynman diagrams contributing to the 2— and the 4—point function. We 
start with the correction to the quadratic action given by 


where 


AS'quad[<?^] — MI 2 >co 


(3.7) 


= uo Afi'^{pi,ki)+ vo Afj.‘^{pi,ki)coski Api. (3.8) 

J ki J h\ 

The correction Ap^ipi^ki) is given by the sum of all bubble graphs shown on figure ([6]) 
with an effective vertex given by a combination of the planar vertex —uq and the non 
planar vertex — uocosp A k where p and k are the momenta flowing into the vertex as 
shown on hgure ([7]). The result takes the form 


= ^P% + ^plip{Pi)- 
1 


5pp = uq 


^PNpiPi) = Vo 


Iki kj 
Iki kl 


k‘l + p? + 5p?p + 5p\ip{ki) 


1 


cos kl Api- 


(3.9) 

(3.10) 

(3.11) 


+ + 5p?p + 6 plfp{ki) 

Now we discuss the 4—point function. The full correction to the quartic action in the 
large N limit is given by 

ASinM = -^ < Ml >co 

= / ••• / ^a{pi)^aiP2)4>b{P3)MPi) AV {pi, P2, Ps, Pi) ■ (3.12) 

Jpi Jp4 

The leading Feynman diagrams in the large N limit contributing to the correction AV{pi, .■,Pi) 
are shown on figure Q. For simplicity, we skip writing them explicitly. 

The full action: the classical+the complete quantum corrections in the large N limit 
for the quadratic and quartic terms is therefore given by 


-S'eff = / Mp){p‘^ + pIs)M-p) 

Jv<oA 


/p<pA 

+ ... / ^a{Pl)MP2)MP3)MPi) i2TT)'^6‘^{pi + ...+P 4 ) Ves{Pl,P2,P3,Pi)- 

J pi<pA J p4<pA 

(3.13) 


The definition of A/x^g and t4fr are obvious. 
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3.2 The Noncommutative Wilson-Fisher Fixed Point 


Next, we apply the Wilson renormalization group recursion formula to get an explicit 
expression of this action. We will assume for simplicity that uq = vq. After some calcula¬ 
tion we obtain 


5, 


quad 



+ Ap?{p) 


(f>ai-Ap). 


(3.14) 


5'int = / ... / ^a{Api)^a{Ap2)MAp3)^b{Ap4) {2Tr)'^6‘^{pi + ...+P4:) Ves{Api,..,Ap4). 

7pi<p P4^P 

(3.15) 


/i^, uo and 9 are dimensionless parameters, viz = p?/A?, uq = uq/A^ and 6 = 6 A‘^. 
The corrections Ap"^ and are now given by (with Vd = l/A'^) 


Ap^(p) = 2-^VdC2(a). 
I + p^ 


(3.16) 


14 fr(Api,..,Ap 4 ) = 1^1 - 4uol/rfA^(A)C4(a)^ cosA^ ^^ cosA^ ^^ 

, 1 .2P1 AP 3 +P 2 AP 4 , .2P1 AP4+P2 AP 3 

+ - cos A-h cos A - 

2 V 2 2 


.(3.17) 


C 2 and C 4 are the connected two-point and four-point functions of the zero-dimensional 
vector model which is given by the functions 


C2(a) = — 1 = i_a + 2a^-5a^ + 


2 a 


(3.18) 


^ / \ 1 — (1 4fl) 2 93 

(94 a =-^- - — = l-3a + lOa^ - 350^ -b ... 

2 a 


The effective coupling a is defined by 
VdUQ 


a = 


( 1 +/ 22)2 


1 + Xd-i{e) ] , Xd-i{x) = 


(3.19) 




The renormalization group equations which follow from the above effective action are 
given by the equations 


-p^ = p-^ -p^ + t{l +-p^) 


- 2 \ (1 + 4 a) 2 — 1 


(3.21) 


^0 + ^0 = 2/0 uo[l — t 


9 — (1 + 4a) 2 


(3.22) 
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Uq-Vq = - 


p ‘^uot{l - (1 + 4a) 


(3.23) 


We have set 


2 

l + lrf_i(0) 


P 9^ 

^ 8d'^{d + 2) ^ ■■■■ 


(3.24) 


We will also use the notation 

^_^ (l + 4a)i -1 
2 


(3.25) 


A non-gaussian fixed point is given hy = p? = pi and Uq + Ug = 2uo = 2u* or 
equivalently 




2 

* 


T, 

- 1 - T* • 


(3.26) 


= 1 — t 


1 — (1 + 4a*) 2 


(3.27) 


The fixed point near e ~ 0, for any value of the dilatation parameter p, is given by 

elnp 2 elnp elnp 

a* = -— , = —2 - 7 ’ 

t — 1 ItVd 

For 0 < e < 2, the critical values of a and T are given by 


(1 - z){t -l + z) 

= (t_2 + 2.p 


(3.29) 


= t 


1 -P^ 

t - 2 + 2 p^ ■ 


(3.30) 


We must clearly have p^ > 1 — t/2 and < 1 for T* to be positive definite. We get 
therefore 


,__ 

2(1 -t) + {t- 2 )(p 2 + p^) + 2 p ^+2 ’ 

t (l-p 2 ) 2 (l_^.)(t_l + p.) 

2Frf (2(1 - t) + (t - 2) (p2 + p^) + 2p^+2) 2 ■ 


We can check that we have always a* > 0 and pi < 0. 

The perturbative solution (13.2711 works on the perturbative sheet 1—t/2 < p^ < 1. 
However, there should be no difference between the regions 1 —t/2<p'^<l and 0 < p*^ < 
1 — 1/2 and thus one must analytically continue the above solution to the non-perturbative 
sheet 0 < p^ < 1 — t/ 2 . 

On the perturbative sheet the function f{z) starts from / = 0 at z = 1 then increases 
to oo as z decreases to z = 1 — t/ 2 , whereas on the non-perturbative sheet it starts from 
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/ = {t — — 2)^ > 0 at 2 ; = 0 then increases to 00 as z increases to z = 1 — t/2. The 

function /(z) vansihes at z = 1 — t < 0. We have then 


(3.32) 

The plus sign corresponds to the above perturbative solution. The minus sign in (|3.32|) 
leads to 


2-t-p^ 


I _ ~ ^ 

2 


(3.33) 


By setting z = 2—t—on the second sheet, equation (j3.33|) becomes equation (j3.27p . This 
leads to a* = /(z) where /(z) is given by equation (j3.29p . The first sheet corresponds to 
the interval z E [1 — t/2,1] whereas the second sheet corresponds to z E [0, 1 — t/2]. When 
we continue the solution to the second sheet we observe that the critical coupling constant 
a* does not return to 0 when we take the limit —)• 0. Indeed, /(O) = (t—l)/(t—2)^ > 0. 
As long as 9 is sufficiently small we have t near 1 and as a consequence /(O) is small and 
we get the commutative result. The critical value of T on the second sheet is 


^ t(t- 1 + p^) 

* t-2 + 2p^ ■ 

As a consequence the critical values pi and u* on the second sheet are given by 


(3.34) 


t{t - 1 + p*^) 

~ 2 - 2t + + (i _ 2)(p2 + pE) 2p2+^ ’ 

t (t-l + p^)(l-p^)(l-p2)2 

(2 - 2t + t2 + (t _ 2)(p2 + 2p2+^)2 ■ 


(3.35) 


This is the noncommutative Wilson-Fisher fixed point. 

We observe that for a fixed t the limit of the perturbative fixed point (13.311) when 
^ 1 is (with XjVd = 2'^-i7r‘^/2r(d/2)) 


e In p 2 £ f 

a, = — , p, = , u. = ^. 


(3.36) 


The limit for a fixed t of the non-commutative Wilson-Fisher fixed point (|3.35p when 
p*^ —;► 0 is 

4 


* “ 1 n , P' 
O* — T--rT7T[l + 


{t - 2)2 


t-l 


t-2 


+ -J , 


pI = 


p' 


t - 2 


{t-iy + v t-l (f - 1)2 +1 


P^ + -] , 


u. = T ,„ [1 + ^ + ...]. (3,37) 


2ri((t-l)2 + l)2‘ 't-l 


In contrast with the commutative theory and with the perturbative fixed point, the non¬ 
commutative Wilson-Fisher fixed point is not vanishingly small in the limit p*^ —)• 0 and 
becomes significantly more important as we increase t from 1 to 2, i.e. as we increase the 
noncommutativity 0 from 0 to tt. Indeed, we see that a* —00 when t —)• 2, i.e. when 
we have only one sheet [0,1]. Putting it differently, in the limit t —>■ 2 the two-sheeted 
structure of a* disappears and we end up only with the perturbative fixed point (I3.3ip . 
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Figure 1: Feynman digrams contribnting to the renormalization of the mass parameter and the 
harmonic oscillator conpling constant. 



Figure 2: Feynman digrams contribnting to the wave fnnction renormalization (linear term in 
and also to the renormalization of the mass parameter and the harmonic oscillator coupling 
constant (the = 0 term). 



Figure 3: The eigenvalue A 3 as a function of the dilatation parameter p. 
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Figure 4: The eigenvalues A* as a function of the dilatation parameter p. 



Figure 5: The anomalous dimension rj as a function of the dilatation parameter p. 



—NV(ki, ^ 2 , /ca, kp — NV(k, k, 1,1) = —ug — vgcos(k A 1) 


Figure 7: The vertex of non-commutative 0{N) sigma model. 
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Figure 8: The dressed vertex of non-commutative 0{N) sigma model. 
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